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The electric charge density in the vortex lattice of superconductors is studied within the Ginzburg- 
Landau theory. We show that the electrostatic potential ip is proportional to the GL function, 
(f oc — IV'ool^- Numerical results for the triangular vortex lattice are presented. 



Abrikosov vortices in type-II superconductor consist 
of magnetic flux encircled by a rotating condensate of 
Cooper pairs. The inertial and Lorentz forces both act in 
the centrifugal direction and thus cause a small depletion 
of the charges in the vortex core. The electrostatic po- 
tential generated by this depleted charge supplies a cen- 
tripetal force which balances the centrifugal force . 

While the depleted charge occupies a region compara- 
ble with the vortex core, i.e. a region extending over the 
Ginzburg-Landau (GL) coherence length ^, the screen- 
ing charge extends over a region extending over half the 
London penetration depth A. This charge distribution 
results in an electrostatic field which may be observed 
either by nuclear quadrupole resonance js) or by surface 
scanning as proposed recently by Blatter et al. Q or by 
the Torricelli-BernouUi effect Besides, the charge of 
the vortex core is expected to contribute to the forces 
acting on the vortex. For example, Feigel'man et al. ||] 
speculate that this charge can explain the still puzzling 
sign reversal of the Hall voltage. 

The charge in the asymptotic region away from the 
vortex core has been studied by LeBlanc In this 
region one can benefit from the simplicity of the London 
theory. Within the London theory, the local balance of 
forces acting on the Cooper pair may be written in form 
of the Newton equation. 
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Here the mass m* = 2m and charge e* = 2e are twice 
the values of a single electron. Unlike in common applica- 
tions of the Newton equation (|^) here the velocity is given 
by the London relation, mv = — eA, and Eq. (|^) thus 
specifies the electric field — which maintains this mo- 
tion. Using the identity V = (vV)v = (v^)— vxVxv 
and the London approximation eB — — mV x v, one can 
show that Eq. (|^) is solved by the Bernoulli potential. 
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The quadratic dependence on the velocity shows that the 
electric potential ip decreases over a characteristic length 
A/2. The charge density p — —eV'^'p behaves similarly. 

Obviously, the Bernoulli potential (j^) is not suited for 
the core region, where the local London approximation 



fails. Moreover, close to the core the velocity diverges, 



" ' 27™rr' "^tiere $0 is the quantum of flux and r is the 
distance from the vortex axis. According to Eq. (||) and 
p = —eV'^ip, the total screening charge follows this diver- 
gence, drp oc r~^. The charge of the core estimated 
from the London expression thus crucially depends 
on the estimate of the core diameter and the profile of 
the charge is incorrect in the core region. 

In the core region, a non-local approach like the 
Ginzburg-Landau theory is required |^ . Like the London 
approach, the Ginzburg-Landau theory does not include 
the electrostatic potential explicitly. Here we show how 
this potential can be determined from its effect on the 
motion of Cooper pairs. The Bernoulli potential then 
has to be replaced by the potential 
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with (3 and Voo having the standard meaning in the 
Ginzburg-Landau theory 

The proof of (||) parallels the proof of the Bernoulli 
potential. The Cooper pairs move according to the 
Schrodingcr equation, 
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^ (iW + e*A)^ V + eVV' = , 



(4) 



which replaces the force balance equation (|l|). 

The electrostatic potential ip also causes a mismatch 
between the local Fermi level and the global chemical 
potential that results in a local suppression of Cooper 
pairs. Briefiy, the charge densities of superconducting 
and normal electrons have a ratio different from the un- 
perturbed state, Ps : Pn 7^ pT ■ P'n ^ where the super- 
script 00 denotes values far from the perturbation. Due 
to this mismatch, the share of normal electrons is locally 
increased. In the linear approximation one has 



Pn =p- -—f. 



(5) 



Taking the gradient, one can read as balance between 
the electric force and a kind of osmotic pressure of normal 
electrons, eE = DV{p^/e). In this sense, equation (|^) 
describes the response of normal electrons to the electric 
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field, importance of which has already been noticed by 
London The "diflFusion coefficient" D is identified 
below. 

Finally, the electrostatic potential obeys the Poisson 
equation. 



wave function. With eV^iyS w the Poisson equation (|^) 
simplifies to Eq. (||). 

Now we compare the above result with alternative ap- 
proaches. First, far from the core one can make the local 
approximation of the Ginzburg-Landau equation (IIG) 



eV = p = + pn + 
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With the wave function normalized to the charge density 
of the condensate, ps — e*|V'P, we can use the condition 
of neutrality, e*|V'ooP + piattice = 0, to efiminate the 
normal density from the Poisson equation. 



£VV = e*(|^AP~lV'ooP)--^^. 
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The set of equations and is closed by the Maxwell 
equation, V x V x A = /ij, with the superconducting 
current j = -^Re [^(ihV + e*A)V']. 

Let us suppose that the system is nearly neutral so that 
the left hand side of (0) may be disregarded. In this case, 
the potential can be determined from the suppression of 
the superconducting density. 
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Substitution of (H) into the Schrodinder equation (^) 
yields a closed equation for ip, 

(i?iV + e*A)^V + 4i?|V'PV'-4i?|V'ooPV' = 0. (9) 
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One can see that for AD = /3 and the standard GL rela- 
tion, a = — /?|'0cx)P, Eq. (ph turns into the GL equation. 
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With D = (3/4:, Eq. (|) coincides with Eq. (|). 

Now we show that the charge neutrality is a good ap- 
proximation in the general case. A simple estimate shows 
that the deviation from charge neutrality, p = — eV^tp, 
is much smaller than the increase of the normal electron 
charge density — 4e^/3~^(p. Indeed, the potential varies on 
the scale of the coherence length, V'^ip ~ £,~'^'P, therefore 
the ratio of these two charge densities behaves as 
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Here we have used the BCS value of /3 see Eq. (13), 
and introduced the GL parameter k — A/f. Even for 
extreme values of the high-Tc materials [k ^ 100, Tc ^ 
90 K, the Fermi energy ep ~ 0.08 eV, m ^ 4me, £ — 4£o, 
~ — piattice] the estimate ( |lT|) gives about 10^^. For 
conventional superconductors the ratio ([ll] ) is typically 
about 10~^. Accordingly, the deviation from the charge 
neutrality may be disregarded when one solves for the 
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where the velocity is given by the current, j — e*v|^p. 
One can see that in this limit the electrostatic potential 
found from ([l^) and (||) equals the Bernoulli potential 
(^. We note that the incompressible superconducting 
liquid discussed by London |l| corresponds to the limit 
K ^ oo of the GL theory |M. In this limit, the local 
approximation is exact. 

Second, formula (|^) is similar to the formula proposed 
by Khomskii and Kusmartsev ]lO| , pT[ , (/Jkk = A^/ {Aeep). 
It turns out, however, that from the BCS-Gor'kov rela- 
tions 101, 
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one obtains 
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i.e., our value is three times larger than tpKK- The sim- 
ilarity of ( p^ with the potential of Khomskii and Kus- 
martsev is, however, accidental. While ( p^ ) follows from 
the intertial and Lorentz forces, the potential ipKK results 
from the condensation energy of the pairing. We have ne- 
glected the condensation energy in our treatment. 

In the rest of this paper we show actual charge profiles 
of an isolated vortex and of the vortex lattice. We use 
known solutions of the Ginzburg-Landau equation and 
evaluate the electrostatic potential if from Eq. (^). The 
deviation of the charge density from neutrality then fol- 
lows from the Poisson equation, p = — eV^(p. 

We start the discussion using the isolated vortex model 
suggested by Clem 
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where ~ ^ is given by the minimum of the free energy 
as = Cn^VI - ^o(^v/A)/i^i'(ev/A), where Kq and 
Ki are modified Bessell functions. This wave function 
fails at long distances, r ~ A, where the correct behavior 
follows the local approximation, IV'P = IV'ooP(l + "^2a )' 
but it provides a good approximation in the vicinity of 
the core, r ~ ^, where the dominant charge contribution 
is located. The electric field, E ~ —'^if, obtained from 
(|l5| ) and (|^) is directed radially. 
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The charge density corresponding to (pi 



-1.3 mV, is identical to the maximum poten- 
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can be compared with the re- 
who found a more compHcated 
We note that Blatter et al. 



is depleted in the core, r < ^v, and enhanced outside the 
core, r > i^v Here, as = is the Bohr radius in the 

crystal. 

The second form of i 
suit of Blatter et al. 
factor, pbi ~ p x tt"^- 
approached the problem from the microscopic side using 
the BCS relation between the total charge and the chem- 
ical potential. This charge is moreover screened in the 
spirit of Jakeman and Pike [p^ within the Thomas- Fermi 
approximation. Since ~ ln(?iijJ_D/Tc) ^ 1 to 10, 

their numerical results are close to ours. 

The charge of the vortex core per unit length, Q = 



2n /q" rpdr — 2neE{S,v)^v, results as 
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In a single pancake (a vortex of length d = 1.17 • 10 ^ m) 
in YBazCuaOy-s (YBCO) at 77 K, the depleted charge 
is Q*d = —0.004 e. We have used ob = 0.05 nm and ^ — 
1.91 nm, where the coherence lenght follows from ^ = A/k 
with A = 191 nm and k = 100. With the quasiparticle 
screening in the spirit of van Vijfeijken and Staas p4| , 
this value will be reduced by ps/p = 2\tpoo\'^/N « 1 — 
{T/Tc)^ ^ 0.5 for T 77 K. 

Note that the vortex charge ( |l8| ) has a value inde- 
pendent of the core diameter A similar value, Q — 
— 1.28 X aBe/(4<^^), results for a different model function, 
= IV'ool tanh(r/^), recommended in This seems 
to contradict the expectation from the local approxima- 
tion, where the core diameter plays a dominant role. In 
the local approach, the velocity close to the core reads 
V ~ 2mnr ' ^rom the Bernoulli potential (^) and the Pois- 
son equation one finds that inside the cylinder of radius 
r the local approximation predicts a charge 
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With the core radius equal to the coherence length, 
r ^ ^ = h/ ^/— 2m*a, the classical estimate ( [1^ ) repro- 
duces the nonlocal result ( p^ . Although the exact agree- 
ment of Eq. ( p^ with Eq. ( |l8| ) is rather accidental, it 
shows that classical estimates used in the literature are 
reasonable. 

An example of the electrostatic potential for the tri- 
angular vortex lattice is presented in Fig. |l|. One can 
see that in the core region the potential has a shape 
similar to the Clem model. The maximum amplitude. 



tial obtained from the Clem model since the GL function 
IV'P reaches zero at the vortex center in both cases. The 
GL function used in Fig. |l| has been evaluated by the 
iteration method described in p^ . 
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FIG. 1. Electrostatic potential in a superconductor with 
triangular vortex lattice of spacing a. The shape of the 
potential follows the magnitude of the GL function, see 
Eq. (^), which depends exclusively on the reduced induction 
b = B/Bc2 = 0.05 and the GL parameter k = 100. The pre- 
sented magnitude of ip results for A = 191 nm and m = 4mo, 
which corresponds to YBCO at temperature 77 K. 
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FIG. 2. Charge density p = —eV^ifi in a superconductor 
with a triangular vortex lattice. All the parameters are the 
same as in Fig. 0. The charge density is plotted in units 
10~^ electrons per unit cell, with the volume of the unit cell 
0.173 nm"^ corresponding to YBCO. 

The density of charge corresponding to the potential in 
Fig. |l| is shown in Fig. ||. In the core, the charge is nega- 
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tive, i.e., the density of charge carriers is depleted. Out- 
side the core, the charge is positive, i.e. there is charge 
built up which screens the charge of the core. The radius 
of the vortex core (radius at which the charge changes 
sign) found in Fig. ^ is 1.53^, comparable to the value 
1.41^ obtained from the Clem model. The maximum 
of the screening charge appears at 2.4^, while the Clem 
model yields 2^. The maximum of the depleted charge in 
the core centre is —5 • 10~^ electron/cell while the Clem 
model gives — 6.8T0~^ electron/cell. The core charge per 
unit lenght is thus reduced to 85% of the Clem model. 
These numerical results thus show that the approxima- 
tion based on the Clem model is reasonable inside and in 
the vicinity of the core. 

In conclusion, we have shown that the GL function 
can be used to evaluate the electric field in superconduc- 
tors, see Eq. (^, in spite of the fact that the electro- 
static potential is absent in the GL theory. For conve- 
nient estimates of the core charge we have derived the 
simple formula (|l8|). Our theory neglects contributions 
from the condensation energy and the screening by nor- 
mal particles; more work on this is under way. Due to 
the neglected screening, the theory applies only for in- 
termediate temperatures where the GL theory still gives 
reasonable results while the share of the superconducting 
electrons is already of the order of unity. Numerical re- 
sults relevant to YBCO at the temperature 77 K and at 
low magnetic field with well separated vortex cores have 
been presented. 
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